MethylExtract: High-Quality methylation maps and SNV calling from whole genome bisulfite sequencing data

[v2; ref status: indexed, http://f1000r.es/301]Whole genome methylation profiling at a single cytosine resolution is now feasible due to the advent of high-throughput sequencing techniques together with bisulfite treatment of the DNA. To obtain the methylation value of each individual cytosine, the bisulfite-treated sequence reads are first aligned to a reference genome, and then the profiling of the methylation levels is done from the alignments. A huge effort has been made to quickly and correctly align the reads and many different algorithms and programs to do this have been created. However, the second step is just as crucial and non-trivial, but much less attention has been paid to the final inference of the methylation states. Important error sources do exist, such as sequencing errors, bisulfite failure, clonal reads, and single nucleotide variants. We developed MethylExtract, a user friendly tool to: i) generate high quality, whole genome methylation maps and ii) detect sequence variation within the same sample preparation. The program is implemented into a single script and takes into account all major error sources. MethylExtract detects variation (SNVs – Single Nucleotide Variants) in a similar way to VarScan, a very sensitive method extensively used in SNV and genotype calling based on non-bisulfite-treated reads. The usefulness of MethylExtract is shown by means of extensive benchmarking based on artificial bisulfite-treated reads and a comparison to a recently published method, called Bis-SNP. MethylExtract is able to detect SNVs within High-Throughput Sequencing experiments of bisulfite treated DNA at the same time as it generates high quality methylation maps. This simultaneous detection of DNA methylation and sequence variation is crucial for many downstream analyses, for example when deciphering the impact of SNVs on differential methylation. An exclusive feature of MethylExtract, in comparison with existing software, is the possibility to assess the bisulfite failure in a statistical way. The source code, tutorial and artificial bisulfite datasets are available at http://bioinfo2.ugr.es/MethylExtract/ and http://sourceforge.net/projects/methylextract/, and also permanently accessible from 10.5281/zenodo.7144.This work was supported by the Spanish Government [BIO2008-01353 to JLO and BIO2010-20219 to MH], and Basque country 'AE' grant (GB)

The Toll-Like Receptor Gene Family Is Integrated into Human DNA Damage and p53 Networks

In recent years the functions that the p53 tumor suppressor plays in human biology have been greatly extended beyond “guardian of the genome.” Our studies of promoter response element sequences targeted by the p53 master regulatory transcription factor suggest a general role for this DNA damage and stress-responsive regulator in the control of human Toll-like receptor (TLR) gene expression. The TLR gene family mediates innate immunity to a wide variety of pathogenic threats through recognition of conserved pathogen-associated molecular motifs. Using primary human immune cells, we have examined expression of the entire TLR gene family following exposure to anti-cancer agents that induce the p53 network. Expression of all TLR genes, TLR1 to TLR10, in blood lymphocytes and alveolar macrophages from healthy volunteers can be induced by DNA metabolic stressors. However, there is considerable inter-individual variability. Most of the TLR genes respond to p53 via canonical as well as noncanonical promoter binding sites. Importantly, the integration of the TLR gene family into the p53 network is unique to primates, a recurrent theme raised for other gene families in our previous studies. Furthermore, a polymorphism in a TLR8 response element provides the first human example of a p53 target sequence specifically responsible for endogenous gene induction. These findings—demonstrating that the human innate immune system, including downstream induction of cytokines, can be modulated by DNA metabolic stress—have many implications for health and disease, as well as for understanding the evolution of damage and p53 responsive networks

Probing the Functional Impact of Sequence Variation on p53-DNA Interactions Using a Novel Microsphere Assay for Protein-DNA Binding with Human Cell Extracts

The p53 tumor suppressor regulates its target genes through sequence-specific binding to DNA response elements (REs). Although numerous p53 REs are established, the thousands more identified by bioinformatics are not easily subjected to comparative functional evaluation. To examine the relationship between RE sequence variation—including polymorphisms—and p53 binding, we have developed a multiplex format microsphere assay of protein-DNA binding (MAPD) for p53 in nuclear extracts. Using MAPD we measured sequence-specific p53 binding of doxorubicin-activated or transiently expressed p53 to REs from established p53 target genes and p53 consensus REs. To assess the sensitivity and scalability of the assay, we tested 16 variants of the p21 target sequence and a 62-multiplex set of single nucleotide (nt) variants of the p53 consensus sequence and found many changes in p53 binding that are not captured by current computational binding models. A group of eight single nucleotide polymorphisms (SNPs) was examined and binding profiles closely matched transactivation capability tested in luciferase constructs. The in vitro binding characteristics of p53 in nuclear extracts recapitulated the cellular in vivo transactivation capabilities for eight well-established human REs measured by luciferase assay. Using a set of 26 bona fide REs, we observed distinct binding patterns characteristic of transiently expressed wild type and mutant p53s. This microsphere assay system utilizes biologically meaningful cell extracts in a multiplexed, quantitative, in vitro format that provides a powerful experimental tool for elucidating the functional impact of sequence polymorphism and protein variation on protein/DNA binding in transcriptional networks

p53 Transactivation and the Impact of Mutations, Cofactors and Small Molecules Using a Simplified Yeast-Based Screening System

The p53 tumor suppressor, which is altered in most cancers, is a sequence-specific transcription factor that is able to modulate the expression of many target genes and influence a variety of cellular pathways. Inactivation of the p53 pathway in cancer frequently occurs through the expression of mutant p53 protein. In tumors that retain wild type p53, the pathway can be altered by upstream modulators, particularly the p53 negative regulators MDM2 and MDM4. promoter, ii) single copy, chromosomally located p53-responsive and control luminescence reporters, iii) enhanced chemical uptake using modified ABC-transporters, iv) small-volume formats for treatment and dual-luciferase assays, and v) opportunities to co-express p53 with other cofactor proteins. This robust system can distinguish different levels of expression of WT and mutant p53 as well as interactions with MDM2 or 53BP1.We found that the small molecules Nutlin and RITA could both relieve the MDM2-dependent inhibition of WT p53 transactivation function, while only RITA could impact p53/53BP1 functional interactions. PRIMA-1 was ineffective in modifying the transactivation capacity of WT p53 and missense p53 mutations. This dual-luciferase assay can, therefore, provide a high-throughput assessment tool for investigating a matrix of factors that can influence the p53 network, including the effectiveness of newly developed small molecules, on WT and tumor-associated p53 mutants as well as interacting proteins

Effects of Ploidy and Recombination on Evolution of Robustness in a Model of the Segment Polarity Network

Many genetic networks are astonishingly robust to quantitative variation, allowing these networks to continue functioning in the face of mutation and environmental perturbation. However, the evolution of such robustness remains poorly understood for real genetic networks. Here we explore whether and how ploidy and recombination affect the evolution of robustness in a detailed computational model of the segment polarity network. We introduce a novel computational method that predicts the quantitative values of biochemical parameters from bit sequences representing genotype, allowing our model to bridge genotype to phenotype. Using this, we simulate 2,000 generations of evolution in a population of individuals under stabilizing and truncation selection, selecting for individuals that could sharpen the initial pattern of engrailed and wingless expression. Robustness was measured by simulating a mutation in the network and measuring the effect on the engrailed and wingless patterns; higher robustness corresponded to insensitivity of this pattern to perturbation. We compared robustness in diploid and haploid populations, with either asexual or sexual reproduction. In all cases, robustness increased, and the greatest increase was in diploid sexual populations; diploidy and sex synergized to evolve greater robustness than either acting alone. Diploidy conferred increased robustness by allowing most deleterious mutations to be rescued by a working allele. Sex (recombination) conferred a robustness advantage through “survival of the compatible”: those alleles that can work with a wide variety of genetically diverse partners persist, and this selects for robust alleles

Mobile DNA elements in T4 and related phages

Mobile genetic elements are common inhabitants of virtually every genome where they can exert profound influences on genome structure and function in addition to promoting their own spread within and between genomes. Phage T4 and related phage have long served as a model system for understanding the molecular mechanisms by which a certain class of mobile DNA, homing endonucleases, promote their spread. Homing endonucleases are site-specific DNA endonucleases that initiate mobility by introducing double-strand breaks at defined positions in genomes lacking the endonuclease gene, stimulating repair and recombination pathways that mobilize the endonuclease coding region. In phage T4, homing endonucleases were first discovered as encoded within the self-splicing td, nrdB and nrdD introns of T4. Genomic data has revealed that homing endonucleases are extremely widespread in T-even-like phage, as evidenced by the astounding fact that ~11% of the T4 genome encodes homing endonuclease genes, with most of them located outside of self-splicing introns. Detailed studies of the mobile td intron and its encoded endonuclease, I-TevI, have laid the foundation for genetic, biochemical and structural aspects that regulate the mobility process, and more recently have provided insights into regulation of homing endonuclease function. Here, we summarize the current state of knowledge regarding T4-encoded homing endonucleases, with particular emphasis on the td/I-TevI model system. We also discuss recent progress in the biology of free-standing endonucleases, and present areas of future research for this fascinating class of mobile genetic elements

Noncanonical DNA Motifs as Transactivation Targets by Wild Type and Mutant p53

Sequence-specific binding by the human p53 master regulator is critical to its tumor suppressor activity in response to environmental stresses. p53 binds as a tetramer to two decameric half-sites separated by 0–13 nucleotides (nt), originally defined by the consensus RRRCWWGYYY (n = 0–13) RRRCWWGYYY. To better understand the role of sequence, organization, and level of p53 on transactivation at target response elements (REs) by wild type (WT) and mutant p53, we deconstructed the functional p53 canonical consensus sequence using budding yeast and human cell systems. Contrary to early reports on binding in vitro, small increases in distance between decamer half-sites greatly reduces p53 transactivation, as demonstrated for the natural TIGER RE. This was confirmed with human cell extracts using a newly developed, semi–in vitro microsphere binding assay. These results contrast with the synergistic increase in transactivation from a pair of weak, full-site REs in the MDM2 promoter that are separated by an evolutionary conserved 17 bp spacer. Surprisingly, there can be substantial transactivation at noncanonical ½-(a single decamer) and ¾-sites, some of which were originally classified as biologically relevant canonical consensus sequences including PIDD and Apaf-1. p53 family members p63 and p73 yielded similar results. Efficient transactivation from noncanonical elements requires tetrameric p53, and the presence of the carboxy terminal, non-specific DNA binding domain enhanced transactivation from noncanonical sequences. Our findings demonstrate that RE sequence, organization, and level of p53 can strongly impact p53-mediated transactivation, thereby changing the view of what constitutes a functional p53 target. Importantly, inclusion of ½- and ¾-site REs greatly expands the p53 master regulatory network

Community-Driven Metadata Standards for Agricultural Microbiome Research

Accelerating the pace of microbiome science to enhance crop productivity and agroecosystem health will require transdisciplinary studies, comparisons among datasets, and synthetic analyses of research from diverse crop management contexts. However, despite the widespread availability of crop-associated microbiome data, variation in field sampling and laboratory processing methodologies, as well as metadata collection and reporting, significantly constrains the potential for integrative and comparative analyses. Here we discuss the need for agriculture-specific metadata standards for microbiome research, and propose a list of “required” and “desirable” metadata categories and ontologies essential to be included in a future minimum information metadata standards checklist for describing agricultural microbiome studies. We begin by briefly reviewing existing metadata standards relevant to agricultural microbiome research, and describe ongoing efforts to enhance the potential for integration of data across research studies. Our goal is not to delineate a fixed list of metadata requirements. Instead, we hope to advance the field by providing a starting point for discussion, and inspire researchers to adopt standardized procedures for collecting and reporting consistent and well-annotated metadata for agricultural microbiome research